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Abstract. The most puzzling property of the extrasolar planets discovered by 
recent radial velocity surveys is their high orbital eccentricities, which are very 
difficult to explain within our current theoretical paradigm for planet formation. 
Current data reveal that at least 25% of these planets, including some with par- 
ticularly high eccentricities, are orbiting a component of a binary star system. The 
presence of a distant companion can cause significant secular perturbations in the 
orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity 
perturbations can occur. These are known as "Kozai cycles" and their amplitude 
is purely dependent on the relative orbital inclination. Assuming that every planet 
host star also has a (possibly unseen, e.g., substellar) distant companion, with rea- 
sonable distributions of orbital parameters and masses, we determine the resulting 
eccentricity distribution of planets and compare it to observations? We find that 
perturbations from a binary companion always appear to produce an excess of plan- 
ets with both very high (e 0.6) and very low (e SS 0.1) eccentricities. The paucity 
of near-circular orbits in the observed sample implies that at least one additional 
mechanism must be increasing eccentricities. On the other hand, the overproduction 
of very high eccentricities observed in our models could be combined with plausible 
circularization mechanisms (e.g., friction from residual gas) to create more planets 
with intermediate eccentricities (e ~ 0.1 — 0.6). 
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1. Introduction 

As of September 2005, more than 160 extrasolar planets have been 
discovered by radial- velocity surveys 1 . At least ~10% are orbiting a 
component of a wide stellar binary system (Eggenberger et al., 2004). 
In contrast to the planets in our own solar system, one of the most 
remarkable properties of these extrasolar planets is their high orbital 
eccentricities. These high orbital eccentricities are probably not signif- 
icantly affected by observational selection effects (Fischer and Marcy, 
1992). Thus, if we assume that planets initially have circular orbits 
when they are formed in a disk, there must be mechanisms that later in- 
crease the orbital eccentricity. A variety of such mechanisms have been 

For an up-to-date catalog of extrasolar planets, see exoplanets.org or 
www. obspm. f r/encycl/encycl .html. 
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proposed (Tremaine and Zakamska, 2004). Of particular importance 
is the Kozai mechanism, a secular interaction between a planet and a 
wide binary companion in a hierarchical triple system with high relative 
inclination (Kozai, 1962; Holman et al., 1997; Ford et al., 2000). When 
the relative inclination angle io between the orbital planes is greater 
than the critical angle i cr it = 39.2° and the semimajor-axes ratio is suf- 
ficiently large (to be in a small-perturbation regime), long-term, cyclic 
angular momentum exchange occurs between the planet and the distant 
companion, and long-period oscillations of the eccentricity and relative 
inclination ensue. To lowest order, the maximum of the eccentricity 
oscillation (ei jmax ) is given by a simple analytic expression: 



(Innanen et al., 1997; Holman et al., 1997). Note that e max depends 
just on iq. Other orbital parameters, such as masses and semimajor 
axes of the planet and the companion, affect only the period of the 
Kozai cycles. Thus, a binary companion as small as a brown dwarf 
or even another Jupiter-size planet can in principle cause a significant 
eccentricity oscillation of the inner planet. 

Our motivation in this study is to investigate the possible global 
effects of the Kozai mechanism on extrasolar planets, and its potential 
to reproduce the unique distribution of observed eccentricities. In prac- 
tice, we run Monte Carlo simulations of hierarchical triple systems. We 
have tested many different plausible models and broadly explored the 
parameter space of such triple systems. 



The purpose of our study is to simulate the orbits of hierarchical triple 
systems and calculate the probability distribution of final eccentricities 
reached by the planet. For each model, 5000 sample hierarchical triple 
systems are generated, with initial orbital parameters based on various 
empirically and theoretically motivated distributions, described below. 
Our sample systems consist of a solar-type host star, a Jupiter-mass 
planet, and a distant F-, G- or K-type main-sequence dwarf (FGK 
dwarf) or brown dwarf companion. The possibility of another giant 
planet being the distant companion is excluded since it would likely be 
nearly coplanar with the inner planet, leading to very small eccentricity 
perturbations. 

The initial orbital parameters of the triple systems are randomly 
generated using the model distributions described in Table I. In this 




(1) 



2. Methods and Assumptions 
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Table I. 



Parameter 


Model Distribution Function 


Ref. 


T-Tnat atar IV/Taaa ry-i ^ ( A/f^, ^ 
nUBL-bLdl IViaoo y+V-LQ ) 


UI111UI III 111 U.y - l.O IvlQ 




Planet Mass mi (Mj up ) 


uniform in log mi, 0.3 — 10Mj up 


[1] 


Secondary Mass...m2 (Mq) 


^(g = m 2 /mi)~ exp | ~ (g ~° 5 23) | 


[2] 


Semimajor Axis....ai (AU) 


uniform in log 01, 0.1 — 10 AU 


[1], [3] 


of Planet 






Binary Period P 2 (days) 


/(logP 2 )~exp{-^^- 48)2 } 


[2] 


Eccentricity of Planet. ..ei 


10" 5 




Age of the System. ..to 


uniform in 1 — 10 Gyr 


[4] 



[1] Zucker and Mazeh (2002), [2] Duquennoy and Mayor (1991), [3] Ida and Lin 
(2004), [4] Donahue (1998) 



Table II. 



Model 


£12,FGK (AU) 


02.BD 


a (AU) 




e b 2 


BDs c 


A 


using P 2 , < 2000 


100 


- 2000 


10~ 5 


- 0.99 


5% 


B 


using P 2 , < 2000 


100 


-2000 


10~ 5 


- 0.99 


10% 


C 


using P 2 , < 2000 


100 


-2000 


10~ 5 


- 0.99 


20% 


D 


using P 2 , < 2000 


100 


-2000 


10~ 5 


- 0.99 


30% 


E 




10 


- 2000 


0.75 


- 0.99 


100% 


F 




10 


-2000 


0.75 


- 0.99 


5% 



a uniform in logarithm 

b all from thermal distribution, P(e 2 ) = 2e2 
c the fraction of brown dwarfs in 5000 samples 



paper, we present six models, each with different initial conditions that 
are listed in Table II. 

For the calculation of the eccentricity oscillations, we integrated 
the octupole-order secular perturbation equations (OSPE) derived in 
Ford et al. (2000). These equations also include GR precession effects, 
which can suppress Kozai oscillations. As noted by Holman et al. (1997) 
and Ford et al. (2000), when the ratio of the Kozai period (Pkoz), to 
the GR precession period (-Pgr) exceeds unity, the Newtonian secular 
perturbations are suppressed, and the inner planet does not experience 
significant oscillation. 

Figure 1 shows typical eccentricity oscillations in two different triple 
systems. One contains a distant brown dwarf companion and the other 
a solar-mass stellar companion. The two systems have the same initial 
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Figure 1. Eccentricity oscillation of a planet caused by a distant brown dwarf 
companion (M = 0.08M Q , solid line) and by a main-sequence dwarf companion 
(M = 0.9M Q , dotted line). For both cases, the mass of the planet host star 
mo = IMq, the planet mass mi = lMj, the planet semimajor axis a\ — 2.5 AU, 
the semimajor axis of the companion ai = 750 AU, the initial eccentricity of the 
companion = 0.8, and the initial relative inclination io = 75°. Note that ei jmax 
is the same in both cases, as it is dependent only on io, but the smaller mass of a 
brown dwarf companion results in a much longer oscillation period fkoz • 

orbital inclination (io = 75°), and we see clearly that the amplitude of 
the eccentricity oscillation is about the same but with a much longer 
period Pkoz for the lower mass companion. 

To find the final orbital eccentricity distribution, each planetary 
orbit in our systems is integrated up to the assumed age of the sys- 
tem (to), and then the final eccentricity (ef) is recorded. The results 
for representative models are compared to the observed eccentricity 
distribution in §3. For more details, see Takeda and Rasio (2005). 

3. Results for the Eccentricity Distribution 

Figure 2 shows the final eccentricity in various models. Each model 
is compared with the distribution derived from all the observed single 
planets with a\ > 0.1, from the California & Carnegie Planet Search 
Catalogue. The statistics of the final eccentricities for our models and 
for the observed sample are presented in Table 3. Models A-D represent 
planets in hierarchical triple systems with orbital parameters that are 
broadly compatible with current observational data and constraints 
on stellar and substellar binary companions. All the models produce a 
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Figure 2. Final cumulative eccentricity distributions {top) and normalized proba- 
bility distributions in histogram [bottom). Four models with different fractions of 
brown dwarf and stellar companions. Initial inclinations (io) are distributed uni- 
formly in cosio (left). Two extreme models where all the binary companions have 
orbits inclined by more than 40° (right). 

Table III. Statistics of Eccentricity Distributions 



Model Mean First Quartile Median Third Quartile 

Observed 0.319 0.150 0.310 0.430 

A 0.213 0.000 0.087 0.348 

B 0.215 0.000 0.091 0.341 

C 0.201 0.000 0.070 0.322 

D 0.203 0.000 0.066 0.327 

E 0.245 0.000 0.141 0.416 

F 0.341 0.071 0.270 0.559 



large excess of planets with ef < 0.1 (more than 50%), compared to only 
19% in the observed sample - excluding multiple-planet systems. An 
excess of planets which remained in low-eccentricity orbits was evident 
in most of the models we tested. Changing the binary parameters, such 
as separations or frequency of brown dwarf companions, did not change 
this result. 

A major difference between most of the simulated and observed ec- 
centricity distributions in the low-eccentricity regime (e < 0.1) mainly 
arises from a large population of binary companions with low orbital 
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inclination angle. For an isotropic distribution of io, about 23% of the 
systems have io < i cr it, leading to negligible eccentricity evolution. For 
completeness, biased distributions of iq and e2 are tested in model E 
and F, as an attempt to achieve the best possible agreement with the 
observations. With all the binary companions having sufficient incli- 
nation angles, model F shows a better agreement with the observed 
sample in the low-eccentricity regime. However, the number of planets 
remaining in nearly circular orbit (e < 0.1) is still larger than in the 
observed sample. Moreover, model F produces the largest excess of 
planets at very high eccentricities (e > 0.6). Note that these extreme 
models are clearly artificial, and our aim here is merely to quantify how 
large a bias would be needed to match the observations "at any cost." 



4. Summary and Discussion 



In most of our simulations, too many planets remain at very low orbital 
eccentricities. The fraction of planets with e < 0.1 is at least 25% in 
our models, but only ~ 15% in the observed sample. There are several 
reasons for this overabundance of low eccentricities in our models. First, 
the assumption of an isotropic distribution of io automatically implies 
that 23% of the systems have io < *crit> resulting in no Kozai oscillation. 
This fraction already exceeds the observed fraction of nearly circular 
orbits (ei < 0.1) which is ~ 15%. Systems with sufficient initial relative 
inclination angles still need to overcome other hurdles to achieve highly 
eccentric orbits. If many of the binary companions are substellar or in 
very wide orbits, Kozai periods become so long that the eccentricity 
oscillations are either suppressed by GR precession, or not completed 
within the age of the system (or both) . This can result in an additional 
15%-40% of planets remaining in nearly circular orbits. Even when the 
orbits of the planets do undergo eccentricity oscillations, there remains 
still 8-14% that simply happen to be observed at low eccentricities. 
Thus, our results suggest that the observed sample has a remarkably 
small population of planets in nearly circular orbits, and other dynam- 
ical processes must clearly be invoked to perturb their orbits. Among 
the most likely mechanisms is planet-planet scattering in multi-planet 
systems, which can easily perturb eccentricities to modest values in the 
intermediate range ~ 0.2 — 0.6 (Rasio and Ford, 1996; Weidenschilling 
and Marzari, 1996; Marzari and Weidenschilling, 2002). Clear evidence 
that planet-planet scattering must have occurred in the v Andromedae 
system has been presented by Ford, Lystad, & Rasio (2005). Even in 
most of the systems in which only one giant planet has been detected 
so far, the second planet could have been ejected as a result of the 
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scattering, or it could have been retained in a much wider, eccentric 
orbit, making it hard to detect by Doppler spectroscopy. 

In the high-eccentricity region, where e\ k, 0.6, our models show 
much better agreement with the observed distribution. The Kozai mech- 
anism can even produce a small excess of systems at the highest ec- 
centricities (ei > 0.7), although it should be noted that the observed 
eccentricity distribution in this range is not yet well constrained. It 
is evident that the observed planets are rather abundant in interme- 
diate values of eccentricity. The Kozai mechanism tends to populate 
somewhat higher eccentricities, since during the eccentricity oscillation 
planets spend more time around ei jmax than at intermediate values. 
However, this slight excess of highly eccentric orbits could easily be 
eliminated by invoking various circularization processes. For example, 
some residual gas may be present in the system, leading to circulariza- 
tion by gas drag (Adams and Laughlin, 2003). In another scenario, 
decreased periastron distances can consequently remove the orbital 
energy of the planet by tidal dissipation. This mechanism, referred 
to as "Kozai migration", was proposed by Wu and Murray (2003) to 
explain the orbit of HD80606 b. Kozai migration can also circularize 
the planetary orbit. It is worth to note that the only three massive 
hot Jupiters in the observed sample, rBoo b, GJ86 b and HD195019 b 
(Msini > 2Mj up , P < 40 days) are all in wide binary systems (Zucker 
and Mazeh, 2002). Their tight orbits with low eccentricity can be a 
consequence of wider orbits with small periastron distances, initially 
invoked by the Kozai mechanism. 

Clearly, even by stretching our assumptions, it is not possible to ex- 
plain the observed eccentricity distribution of extrasolar planets solely 
by invoking the presence of binary companions, even if these compan- 
ions are largely undetected or unconstrained by observations. However, 
our models suggest that Kozai-type perturbations could play an impor- 
tant role in shaping the eccentricity distribution of extrasolar planets, 
especially at the high end. In addition, they predict what the eccen- 
tricity distribution for planets observed around stars in wide binary 
systems should be. The frequency of planets in binary systems is still 
very uncertain, but the search for new wide binaries among exoplanet 
host stars has been quite successful in the past few years (e.g., Mugrauer 
et al, 2005). 
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